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Abstract 

We investigate possible contributions to future atomic parity nonconservation (PNC) experi¬ 
ments from parity admixtures between single-electron atomic states with total angular momentum 
j = 3/2. We develop new formalism for studying these admixtures between atomic p and d states, 
which enter only when finite nuclear size effects are considered and have been neglected in the lit¬ 
erature to date. We use analytic approximations to provide an order-of-magnitude estimate of the 
contribution from these admixtures, and identify a dimensionless ratio which sets the scale of the 
correction. Using realistic numerical wavefunctions in Ba"*" we confirm the results of our analytic 
expressions, and conclude quantitatively that these novel admixtures are likely to be negligible in 
essentially all cases. 
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INTRODUCTION 


In recent years, precision experiments B, i measuring parity nonconservation (PNC) 
in atomic transitions have developed into low-energy tests of the Standard Model which 
complement high energy collider experiments. PNC effects in atomic transitions arise due 
to neutral-weak interactions between electrons and the nucleus. Atomic PNC experiments 
are sensitive to the radiatively corrected weak charge of the nucleus, which has been 
calculated to one-loop in the Standard Model 0 S 0 and is given at tree level by = — A^- 1 - 
Z(1 — 4sin^6'„). Because radiative corrections enter into low- and high-energy observables 
differently, these precision atomic PNC experiments provide a search for possible physics 
beyond the Standard Model which is complementary to more conventional accelerator based 
searches. 

It is well known B; IH that the experimental observable must be corrected for hnite 
nuclear size when compared to the Standard Model prediction due to variation of the atomic 
electron wavefunctions over the extent of the nucleus. It has been shown HH 

the transition amplitude between initial and final atomic states can be formally factorized 
to separate nuclear and atomic physics contributions. The nuclear structure correction 
to atomic PNC is non-negligible (about 4% in ^^^Cs) and grows as (Za)^, but is reliably 
calculable. Differences between neutron and proton distributions modify this slightly 0 . 0 . 

ITTI. 


The most precise atomic PNC experiment to date involves parity admixtures of initial 
and hnal states with total angular momentum j = 1/2, but there are several other PNC 
measurements which have either been performed p| or are currently in progress [|I^ [13 


looking for PNC effects in transitions involving j = ?)/2 final states. To date, all calculations 
in the literature iiHl, pTSf consider only contributions from parity admixtures between 
S'i /2 and Pi /2 states because the D 3/2 state has vanishing amplitude and derivative at the 
origin. In principle, there could be an additive contribution to the PNC transition amplitude 
arising from opposite parity admixtures to the j = 3/2 hnal state which is induced by a 
hnite sized nucleus. Motivated by the scale of the hnite nuclear size correction io s — p 
mixing amplitudes, we have quantihed the ehect oi p — d mixing when hnite nuclear size 
ehects are taken into account. 

We begin this work by developing formalism necessary for studying parity admixtures 
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between -P 3/2 and -D 3/2 states in future atomic PNC experiments. We define a dimensionless 
ratio which sets the scale of the additive contribution due to p — d mixing relative to the 
dominant s — p mixing term. We then use approximate analytic calculations to provide 
an order-of-magnitude estimate of this ratio. Using realistic numerical wavefunctions, we 
confirm our analytic estimates in Ba^, one of the atomic systems in which PNC effects 
involving a final D-state may be measured |]T2 . 


FORMALISM 


The parity violating electron-nucleon interaction is dominated by exchange of Z°-bosons 
and can be written in terms of vector and axial-vector currents for electrons and nucleons: 


-f^PNC — + VeA 


/e^lN. 


( 1 ) 


The second term depends on the orientation of the nuclear spin and usually amounts to 
at most a few percent of the first, in part because the vector electron-Z° coupling is small 
compared its axial-vector counterpart. Additional suppression arises because the nuclear 
spin-dependent currents do not add coherently like the vector nucleon currents, which are 
independent of the spin of the nucleus. The nuclear spin independent (NSID) piece of 
the interaction is proportional to and can be isolated experimentally by averaging over 
hyperhne states in the transition, so henceforth we neglect the second term in Eq. Q. 

Due to the small effective momentum transfer associated with atomic PNC observables, 
the three-vector part of the nuclear current in the hrst term of Eq. |T] is highly suppressed 
relative to the charge (p = 0) component. Keeping only the dominant charge component, 
we have [|, [T^, 

( 2 ) 


(i|i/pNc,i|j) 


2V2 

Gp 

'2V2 


dVV'ei(r)7^'0ei(r)Pw(r) 


Here Gp is the Fermi constant, is the weak nuclear charge density, Cij{Z) depends on 
the full multi-electron wavefunction, and contains atomic structure effects including many- 
body correlations for a point nucleus calculation; Af = '0ei(O)7^'0ej(O) is a normalization 
factor for the single-electron axial transition matrix element evaluated at the origin; is the 
correction factor for hnite nuclear size. The remaining factor, is the experimentally 
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determined weak neutral charge of the nucleus which is compared directly with the Standard 
Model prediction. The factorization in the second line of Eq. ^ can be used because the PNC 
matrix elements only depend on the electronic wavefunctions over the extent of the nucleus 
where binding energies can be neglected with respect to the Coulomb potential. Therefore, 
apart from an overall normalization factor, all of the axial electronic matrix elements which 
contribute are essentially identical. 

Of the atomic PNC experiments involving j = 3/2 hnal states, the proposed measure¬ 
ment of PNC in a single trapped Ba"'' ion is the theoretically simplest, since it involves 
a single valence electron outside of a tightly closed shell. The others involve at least two 
valence electrons, leading to additional considerations such as conhguration mixing which 
must ultimately be taken into account. In the proposed Ba’*' experiment, atomic PNC effects 
will be measured by studying the allowed E2 transition between the 6 S 1/2 ground state and 
the 5 D 3/2 excited state. This transition has a small additional PNC-induced El amplitude 


given by 112 


= {5Dy2,m'\El\6Sy2,m) 




E 


(5T>3/2, m'\El\nPi/ 2 , m''){nPii 2 , m'^|TrpNc,i|65'i/2, m) 

fheSi/j - ^nPi /2 

( 5 T) 3 / 2 , m'|hfpNc,i|’^-P 3 / 2 , rn") (nP 3 / 2 , m"|Pl| 6 S'i/ 2 , m) 


(3) 


W. 


5D3/2 




3/2 


_ cPNC I cPN 

— ^S—P — 


p—d 5 


where Tg™*" and carry implicit dependence on m\ m. The Aj = 0 selection rule for 

matrix elements of TTpnc,! has already been explicitly used in determining the intermediate 
states present in the sums in this expression BSP- Here we have separated the PNC- 
induced transition amplitude into contributions arising from parity admixtures of the initial 
s-state ^^cl the hnal d-state 

In initial estimates of the size of the PNC-induced amplitude for the Ba’*' and similar 
1 ^ [ 1 ^ only amplitudes arising from s — p mixing were included, with the 


experiments 


assumption that p and d states are not mixed by iPpNC,i- For a point-like nucleus this is an 
exact description of the system (to hrst order in TTpNc,!)- When hnite nuclear size effects 
are included, however, the validity of neglecting contributions from p — d mixing is less clear. 
As discussed in Refs. I, 0, 0, 0, the neutral weak interaction is a contact interaction 
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between nucleons and the transition electron in atomic PNC. Therefore the matrix elements 
of depend on electronic wavefunctions evaluated over the extent of the nucleus. 

Upper- and lower-component Dirac radial wavefunctions vanish as r —0 for j > 3/2, so 
all PNC matrix elements in the second sum in Eq. ^ vanish identically for a point nucleus. 
For a hnite-sized nucleus the hrst term in Eq. ^ still dominates, but it has already been 
shown in Refs. E30 that hnite nuclear size effects can modify at the 5 — 10%-level, 
depending on Z. We now investigate the conditions under which might contribute at 
levels approaching the hnite nuclear size correction to in Eq. p. 

Since is the dominant piece it must be calculated quite accurately, including overall 
normalizations which require sophisticated atomic many-body calculations 01- In order 
to study the contribution from p — d mixing, we factor out this dominant term, leaving a 
small dimensionless correction factor which we estimate roughly: 


'PNC 


where 


:Tg™C [1+Pp_d], 

(4) 

yPNC 

« ga 

) 

III 

1 

a 

(5) 


Formally, this ratio involves sums over transition matrix elements of i^pNcq cind El operators 
between full many-electron atomic states including correlation effects, the calculation of 
which is beyond the scope of the present work. Since our goal is an order-of-magnitude 
estimate for the correction, however, we approximate the transition matrix elements using 
single-electron atomic states of nominally good parity. We expect this to be a reasonable 
approximation for Ba^ because the transition involves a single valence electron outside of a 
tightly bound core. Using the solutions to the Dirac central potential problem with dehnite 


parity [1^ , we perform the angular integration explicitly [|^ , reducing these single-electron 
matrix elements to expressions containing simple radial integrals. We then express the PNC 
matrix elements as 


(?7,Pi/2,m"|iPpNc|6S'i/2,m) = -i 

(5D3/2, m'|PpNc|^-P 3 / 2 , ^m'm" 


t-PNC 

^F --PNC 

2\/2 


( 6 ) 

(7) 
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and the El matrix elements as (taking the ^-component for simplicity) 

n-E\ 

-L't: 


2\/2 

{bD^/2,m'\El^\nPi/2,m") = e h 


m'm" -^5dnp 


2\/2 

{nP3/2,m"\El^\6Si/2,m) = ^ e 6, 


'j~El 

m"m -^np6s ' 


(8) 

(9) 


With these simplihcations applied to Eqs. ^ and onr correction factor arising from p — d 
mixing can be written 


T^p-d — 


lV„p. 


( 10 ) 


In these expressions we hnd it convenient to dehne transition radial integrals, Xnp^n 2 i 2 J 
because phases, numerical constants, and all dependence on magnetic quantum numbers 
cancel in the ratio. 

We now turn to evaluating the transition radial integrals of interest. In terms of single¬ 
electron radial wavefunctions, the integrals arising from the matrix elements of i?pNC,i take 
the form 


XS? = j drp^r) 
[ drp^{r) 


G. 


^Pl/2 (^)-^6 si/2 (^) Epi /2 ij')^dSi/2 (^) 


^5^3/2 {^)Pnp^/2{^) - F^d 

3/2 {^)Gnp^/2{^) 


( 11 ) 

( 12 ) 


These integrals constitute the hnite nuclear size corrections to the s — p and p — d mixing 
contributions to atomic PNC, with the weak charge density of the nucleus (p„ = —Np^ + 
Z(1 — 4sin^6*w)pp) integrated against a radial folding function. These expressions contain 
essentially all dependence of 7Zp-d on nuclear physics. The transition radial integrals arising 
from the El matrix elements are given by H 


iSnp = I * ^ 5 * 


n-El _ _ 

'^npQs 


^/2X)Gnp,/2{r) 
dr Gnp,^^{r)EGs,^^{r). 


(13) 

(14) 


The parity-allowed El transition strengths and energy denominators present in the ratio, 
TZp-d, are generally of similar orders of magnitude |K, 0- Therefore we expect the scale 
of the correction due to p — d mixing to be set by the PNC transition matrix elements. We 
now look for analytic approximations to the transition radial integrals in Eqs. o and 0 
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APPROXIMATE ANALYTIC RESULTS 


Although the integrals in Eqs. 0 and |T^ formally extend over all space, the weak charge 
density only contributes appreciably for r < 10 fm. Therefore the PNC transition radial 
integrals can be estimated reliably by considering the form of the electron wavefunctions in 
the vicinity of the nucleus. In this region the electronic potential is dominated by the nuclear 
Coulomb potential, with screening and correlation effects negligible [^]. Also, the electronic 
binding energies can be safely neglected with respect to the potential in the vicinity of the 
nucleus, as discussed in Refs. 1^. Using power series solutions of the Dirac radial 
equations []^ under these minimal assumptions, we find that the PNC transition radial 
integrals can be expressed as 


jPNC _ 4 

-^npSs 


-N 


1 - g0o(0o - 2me)(r^)n H- 


+Z(1 - 4sin^ 9^) 


>o(0o - 2me)(r2)p H- 


1 

47r 


+ Z{1 - 4 Sitl^ »„) {P)p + ■••}. 


(15) 


(16) 


In obtaining these results we have expanded the realistic nuclear Coulomb potential in a 
Taylor series about the origin, V (r) ~ 0o + 02?"^ + ■ ■ ■. Here 0o is the leading term in 
the expansion of the potential, and Anij is the coefficient of the leading term in a power 
series expansion of the normalized single-electron wavefunction, \nlj,m). We note that the 
quantities in square brackets in Eq. are the finite nuclear size corrections defined in 
I0|, Q. and that the p — d transition radial integrals, vanish for a point 


Refs. 


nucleus as expected. 

Because we are looking for an order-of-magnitude estimate of TZp-d, details like the ~ 5% 
finite nuclear size correction to and possible differences in neutron and proton mean- 

square radii are unimportant at this stage. We further simplify the PNC transition radial 
integrals by setting the finite nuclear size correction equal to unity in Eq. and replacing 
(r^)n and {r‘^)p in Eq. |T^ with the nucleon mean-square radius With these approxi¬ 

mations we find 


7 -PNC 

-^np6s 

7 -PNC 

-^5dnp 


-^^npi/ 2 ^esi /2 + X(1 — 4sin^ 0w)} 


dvr 

1 

dvr 


i5d3/22i„p3/2 {-N + Z{1 - dsin^^w)} (r^)N. 


(17) 


AkJ Ar, 
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It is important to note that the Ani/s in Eqs. and |T^ are dimensionful constants, 
and that the constants for j = 1/2 and j = 3/2 have different units. Both 

arise from matrix elements of the same operator, i^pNcp; so the ratio must be di¬ 
mensionless. This expansion of the integrals for small r shows, therefore, that the ratio 
(^ 5 ^ 3 / 2 ^?up 3 / 2 )/(^fi 2 Pi/ 2 ^ 6 si/ 2 ) should havo dimension where L is a length scale rele¬ 

vant to the normalization of atomic bound state wavefunctions. 

This dependence on a characteristic atomic length scale is indicative of the fundamental 
difference between estimating contributions from p — d mixing and studying the sensitivity 
of the Boulder ^^^Cs experiment to the spatial neutron distribution [^, |l^. In the Boulder 
experiment [|T|, only parity admixtures between s and p states are present, so a formal 
factorization of the transition amplitude (Eq. ^ is possible. This factorization gives separate 
multiplicative factors dependent on contributions from nuclear and atomic physics which can 
be calculated independently. In evaluating 7^p_d no such factorization is possible, so we must 
estimate the normalization constants Anij- Plugging Eqs. 0 into our expression for the p — d 
mixing correction factor, Eq. 0 , we now have 


TZn- 


p-d 


2^np6s/ (^^5D3/2 


(18) 


"^6^1/2 ^np^/^'^idnp! (^6Si/2 “ W^nPi/ 2 ) 

Here, the sums in numerator and denominator should be of similar orders of magnitude |T^, 


17|, so we expect that the overall scale of the correction factor will be given by the ratio out 


in front. We will conhrm this expectation numerically for Ba’*' in the next section. 

In order to estimate the coefficients A^ij in Eqs. ^ and normalized electronic wave- 
functions satisfying the condition {^nij + ^nij) — ^ must be constructed and expanded 
for small r. Hence, calculation of the normalization coefficients at the origin requires knowl¬ 
edge of the electron wavefunctions over atomic distance scales where shielding and correlation 
effects become important. Such calculations are beyond the scope of the present work. We 
note, however, that our correction factor 77.p_d will depend on the ratio of (p^)n to the rel¬ 
evant atomic length scale, and that all atoms are roughly the same size. We therefore hrst 
look to hydrogenic electron wavefunctions for analytic expressions to estimate the Anij and 
study the dependence on the characteristic parameters of the problem. Fractional deviations 
between (upper-component) relativistic and nonrelativistic radial wavefunctions are gener¬ 
ally of order {Za^, except where one of the wavefunctions has a node ||^. We therefore 
look at the normalization of nonrelativistic hydrogenic wavefunctions at the origin for order- 















of-magnitude estimates of the Anij's. The coefficient of the leading term in an expansion of 
the nonrelativistic hydrogenic wavefunction as r —> 0 is given by [| 
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Cnl = 


/ 2Z \ 

{n + l)\ 

1/2 ^ 

\nao) 

to 

1 

1 

(2/ + 1)! 


(19) 


where Oq, the Bohr radius, gives the length scale relevant to atomic wavefunction normal¬ 
izations. Replacing the in Eq. by the nonrelativistic CniA we hnd 


7^p_d ~ 0.04 


Z"(r=)N 


3/2 1 


1 X 10 
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( 20 ) 


E„C'„^fi/(lVes.,.-W'„P,p)’ 

(using Z = 56 and = 4.84 fm for Ba’*'). Based on this rough estimate, it appears 

that the ratio of nuclear to atomic length scales predicted by the analytic approximations 
serves as the primary suppression factor in considering p — d mixing. Thus far we have made 
general arguments that the scale of TZp-d should be set by the ratio of PNC radial integrals 


in Eq. without estimating the El radial integrals or energy denominators. Physically, 
we know that El transition strengths and binding energies for the states of interest here are 
generally within an order of magnitude of one another [K .0- Theoretically, these quantities 
cannot be reliably estimated in the context of our approximations because they depend on 
long range atomic physics. In the next section we use realistic numerical wavefunctions 0 
to evaluate all quantities in Eq. 0 as a check of our approximate analytic results. 


NUMERICAL RESULTS 


The calculations presented in the previous section are instructive in quantifying the mech¬ 
anism of suppression for parity admixtures of atomic states with j > 3/2, but are admittedly 
crude. We now wish to test the validity of the order-of-magnitude estimate as well as investi¬ 
gate whether the energy denominators and El transition radial integrals which we have not 
yet addressed might give rise to enhancements to TZp-d- We have obtained realistic binding 
energies and atomic radial wavefunctions for Ba’*', calculated in a relativistic Dirac-Hartree- 
Fock model I3- Using these tabulated wavefunctions, we have performed the transition 


radial integrals in Eqs. 0 through 0 numerically. The sums in both numerator and de¬ 
nominator of Eq. ITO are dominated by the n = 6 admixed states, but we have included 
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contributions from n = 7, 8 as well. This calculation yields 77p_d ~ 4 x 10“®, in good 
agreement with our crude but analytic estimate. 

In order to study this numerical check on our approximate analytic calculation, we have 
fit the relativistic numerical wavefunctions to polynomials near the origin and extracted the 
coefficients of the leading terms of the upper- and lower-components. Overall wavefunc- 
tion normalizations are not fixed by power series solutions near the origin, but the ratio 
of leading term normalizations of upper- and lower-component wavefunctions is uniquely 
determined by the value of the potential at the origin. Comparing this ratio estimated from 
the htted numerical wavefunctions and from power series solutions for relativistic hydrogenic 
wavefunctions assuming that the total potential at the origin is dominated by the nuclear 
Coulomb potential, we find that the two methods agree to better than 5%. As an additional 
check, we assume that 77p_d is dominated by the n = 6 term and estimate as 

simplihed in Eqs. 0 using both nonrelativistic and htted coefhcients. The nonrelativistic 
hydrogenic estimate gives 1 ■ 10“®, while the estimate from htting numerical wavefunctions 
gives 0.6-10“®. These numerical checks indicate that our crude analytic approximate results 
are in reasonable agreement with a more detailed, full calculation requiring many-electron 
atomic wavefunctions including the effects of correlations and shielding. 


CONCLUSIONS 


These calculations indicate that neglecting possible p — d mixing contributions to atomic 
PNC amplitudes which involve j = 3/2 initial or hnal states is an exceedingly good approx¬ 
imation. We have demonstrated that errors from neglecting p — d mixing effects are well 
below other sub-1% effects which have been neglected to date, including contributions from 
the three-vector nucleon currents, radiative corrections such as one-photon one-Z° exchange 
box diagrams, and parity violating electron-electron interactions. While we have focused 
on the simplest atomic system proposed for experiment, Ba"'', we expect this result to be 
applicable to more complex systems like Ytterbium because the scale of the correction 
is proportional to in all cases. The most likely source of significant enhancement 

of the p — d mixing term would be an exceedingly small energy difference between the the 
admixed p and d states in a heavy atom, which would depend sensitively on the particular 
atomic system. 
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